function [c,ceq] = fminconstr_crs(x,e2,ell,theta,sigma,pi,eta,i,M,A,B)

c = []; % no nonlinear inequality
ceq = sys_crs(x,e2,ell,theta,sigma,pi,eta,i,M,A,B); % the fsolve objective is fmincon constraints